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^ Abstract 
Q 

^ We recap the main features of Minimal Dark Matter (MDM) and 

assess its status in the light of the recent experimental data. The 
theory selects an electroweak 5-plet with hypercharge Y = as a fully 
CIh successful DM candidate, automatically stable against decay and with 

no free parameters: DM is a fermion with a 9.6 TeV mass. The direct 
^ detection cross-section, predicted to be 10~ 44 cm 2 , is within reach 

of next-generation experiments. DM is accompanied by a charged 
> fermion 166 MeV heavier: we discuss how it might manifest. Thanks 

to an electroweak Sommerfeld enhancement of more than 2 orders 



00 

of magnitude, DM annihilations into W + W~ give, in presence of a 
^ modest astrophysical boost factor, an e + flux compatible with the 

PAMELA excess (but not with the ATIC hint for a peak: MDM 
O instead predicts a quasi-power-law spectrum), a p flux concentrated 

at energies above 100 GeV, and to photon fluxes comparable with 
^ present limits, depending on the DM density profile. 



1 Introduction 

The quest for the identification of the missing mass of the universe has been with us since 
many decades now [1]. While explanations in terms of modifications of Newtonian gravity 
or General Relativity become more and more contrived, evidence for the particle nature of 
such Dark Matter (DM) now comes from many astrophysical and cosmological observations. 
Non-baryonic new particles that may fulfill the role of DM have emerged in the latest decades 
within many Beyond the Standard Model (SM) theories, most notably supersymmetry. These 
constructions try to naturally explain the hierarchy between the Electro Weak (EW) scale and 
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the Planck scale and, in doing so, introduce a host of new particles with EW masses and 
interactions. Some of these particles can be good DM candidates (e.g. the lightest neutralino). 

However these approaches to the solution of the DM problem, while still the most popular, 
start facing a sort of 'impasse': (i) The expected new physics at the EW scale has not appeared 
so far at collider experiments: the simplest solutions to the hierarchy problem start needing 
uncomfortably high fine-tunings of the their unknown parameters [2]. (ii) The presence of a 
number of unknown parameters (e.g. all sparticle masses) obscures the phenomenology of the 
DM candidates, (iii) The stability of the DM candidates is usually the result of extra features 
introduced by hand (e.g. i?-parity in supersymmetry) , most often also necessary to recover 
many good properties of the SM that are lost in these extensions (automatic conservation of 
baryon number, lepton number, etc). 

The Minimal Dark Matter (MDM) proposal [3] pursues therefore a different and some- 
what opposite direction: focussing on the Dark Matter problem only, we add to the SM the 
minimal amount of new physics (just one extra EW multiplet X) and search for the minimal 
assignments of its quantum numbers (spin, isospin and hypercharge) that make it a good Dark 
Matter candidate without ruining the positive features of the SM. No ad hoc extra features are 
introduced: the stability of the successful candidates is guaranteed by the SM gauge symmetry 
and by renormalizability. Moreover, due to its minimality, the theory is remarkably predictive: 
no free parameters are present and therefore the phenomenological signatures can be univo- 
cally calculated, e.g. at colliders and for direct /indirect detection, up to the astrophysical and 
cosmological uncertainties. 

The rest of the paper is organized as follows. In Sec. 2 we review the theoretical aspects of 
MDM model. We follow a constructive approach by scanning all the possible choices of quantum 
numbers and selecting the successful DM candidates: insisting on full minimality we will see 
that consistency and phenomenological constraints reject most of the candidates and actually 
individuate only one (the fermionic SU(2) quintuplet with hypercharge Y — 0), of which we will 
study the phenomenology in the following sections. Note that since other candidates can be re- 
allowed by relaxing the request for full minimality (as we discuss in Appendix B), the formulae 
and most of the phenomenological analysis are given in their most general form and can easily be 
adapted to these other candidates. In Sec. 3 we present the computation of the cosmological relic 
density of MDM [4]: imposing the match with the measured density determines the DM mass. 
We also discuss the small mass splitting introduced by loop corrections. In Sec. 4 we review 
the direct detection signatures of the candidate. In Sec. 5 we focus on the indirect detection 
signatures, comparing the predictions ([5]) with the recent data from satellite and balloon 
experiments: this promises to be the most interesting and stringent avenue of confrontation 
with data [6]. In Sec. 6 we review what could be the signatures at colliders. Finally, in Sec. 7 we 
discuss some other possible phenomenological features of the model. In Appendix A we briefly 
review the Sommerfeld enhancement of the annihilation cross section, discussed in more detail 
in [4, 5] . In Appendix B we list the options that open up for model building if the requirement 
of full minimality is relaxed. An executive summary of positive and negative features of the 
model, as well as an outlook towards the implications of future results, is given in Sec. 8. 

2 Construction of the Minimal Dark Matter model 

The MDM model is constructed by simply adding on top of the Standard Model a single 
fermionic or scalar multiplet X charged under the usual SM SUl{2) x Uy(1) electroweak inter- 
actions (that is: a WIMP; it is assumed not to be charged under S77 c (3) strong interactions as 
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Table 1: Book-keeping of the possible Minimal DM candidates and selection of 
successful ones. Quantum numbers are listed in the first 3 columns. The 4th column indi- 
cates some decay modes into SM particles; modes listed in parenthesis correspond to dimension 
5 operators. Candidates with l 7 / are excluded by Direct Detection (DD) searches (un- 
less appropriate non-minimalities are introduced, see App.B), as indicated in the 5th column. 
Candidates with an open decay channel are excluded (unless some other non-minimalities are 
introduced, see App.B), as indicated in the 6th column. Note that, for simplicity, the possibilities 
concerning the SUl{2) 6- and 8-plet are only sketched. Analogously, for the 7-plet we list the 
only interesting candidate. At the end of the game, the fully successful candidates are indicated 
by the shaded background: the fermionic 5-plet with Y = and the scalar 7-plet with Y — 0. As 
for the latter non-minimal scalar quartic couplings are generically present (see Appendix B), 
the former is overall preferred and we will refer to it as the MDM candidate in the paper. 



the bounds are strong on this possibility [7]). Its conjugate X belongs to the same representa- 
tion, so that the theory is vector-like with respect to SUl(2) and anomaly-free. The Lagrangian 
is 'minimal': 

„ _ „ If XiiLp + M)X for fermionic X , , 

j \D^X\ 2 -M 2 \X\ 2 for scalar X {1} 

The gauge- covariant derivative _D M contains the known electroweak gauge couplings to the 
vectors bosons of the SM (Z, W ± and 7) and M is a tree level mass term (the only free 
parameter of the theory). A host of additional term (such as Yukawa couplings with SM fields) 
would in principle be present, but for successful candidates they will be forbidden by gauge 
and Lorentz invariance, as detailed below. 

X is fully determined by the assignments of its quantum numbers under the gauge group: 
the number of its SU(2)l components, n = {2,3,4,5, . . .} and the hypercharge Y. In Table 1 
we list all the potentially successful combinations, as we now proceed to discuss. 
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For a given assignment of n (first column of the table) there are a few choices of the hyper- 
charge Y such that one component of the X multiplet has electric charge Q = T3+Y = (where 
T 3 is the usual 'diagonal' generator of SU(2) L ), as needed for a DM candidate. For instance, 
for the doublet n = 2, since T3 = ±1/2, the only possibility is Y = =Fl/2. For n = 5 one can 
have Y = {0, ±1, ±2}, and so on. We do not consider the case of the n — 1 singlet: lacking 
gauge interactions, even if it is ever produced in the Early Universe it could not annihilate and 
remain with the correct relic amount by means of the standard freeze-out mechanism. 1 
The list of possible candidates has to stop at n < 5 (8) for fermions (scalars) because larger 
multiplets would accelerate the running of the SU(2)l coupling g 2 : demanding that the per- 
turbativity of a^C^O = a^Af) - (6 2 /27r) In E'/M is mantained all the way up to E' ~ Mpi 
(since the Planck scale Mpi is the cutoff scale of the theory) imposes the bound. In this formula 
62 = — 19/6 + c gx{n 2 — l)/36 with c = 1 for fermions, c = 1/4 for scalars and gx is the number 
of degrees of freedom in the multiplet. 

In this list of candidates, those with 7^0 have vector-like interactions with the Z boson 
that produce a tree-level spin-independent elastic cross sections 

r 2 m 2 

cr(DM jV — > DMA/") = c^—^Y 2 (N - (1 - Asi^Z) 2 (2) 

where c = 1 for fermionic DM and c = 4 for scalar DM [8]; Z and N are the number of protons 
and of neutrons in the target nucleus with mass M/v (we are assuming M 3> Mat). This elastic 
cross section is 2 -=-3 orders of magnitude above the present bounds [9, 10] from direct detection 
searches. Unless minimality is abandoned in an appropriate way (that we discuss in Appendix 
B), such MDM candidates are therefore excluded and we will focus in the following on those 
with Y = 0. 

Next we need to inspect which of the remaining candidates are stable against decay into 
SM particles. The fourth column of Table 1 shows some possible decay operators for each case. 
For instance, the fermionic 3-plet with hypercharge Y = would couple through a Yukawa 
operator XLH with a SM lepton doublet L and a Higgs field H and decay in a very short time. 
This is not a viable DM candidate, unless the operator is eliminated by some ad hoc symmetry 
(see again Appendix B). For another instance, the scalar 5-plet with Y = would couple to 
four Higgs fields with a dimension 5 operator XHHH*H* /Mpi, suppressed by one power of the 
Planck scale. Despite the suppression, the resulting typical life-time r ~ Mp\ TeV" 3 is shorter 
than the age of the universe, so that this is not a viable DM candidate. 

Now, the crucial observation is that, given the known SM particle content, the large n multiplets 
cannot couple to SM fields and are therefore automatically stable DM candidates. This is the 
same reason why known massive stable particles (like the proton) are stable: decay modes 
consistent with renormalizability and gauge symmetry do not exist. In other words, for these 
candidates DM stability is explained by an 'accidental symmetry', like proton stability. Among 
the candidates that survived all the previous constraints, only two possibilities then emerge: a 
n = 5 fermion, or a n = 7 scalar. But scalar states may have non-minimal quartic couplings 
with the Higgs field (see Appendix B). We will then set the 7-plet aside and focus on the 
fermionic 5-plet for minimality in the following. 

In summary, the 'Minimal Dark Matter' construction singles out a 

fermionic SU(2)^ 5-plet with hypercharge Y = 
^^We discuss the case of a scalar singlet with non-minimal additional interactions in Appendix B. 
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as providing a fully viable, automatically stable DM particle. It is called 'Minimal DM' (MDM) 
since it is described by the minimal gauge- covariant Lagrangian that one obtains adding the 
minimal amount of new physics to the Standard Model in order to explain the DM problem. 

3 Cosmological relic density and mass determination 

Assuming that DM arises as a thermal relic in the Early Universe, via the standard freeze-out 
process, we can compute the abundance of MDM as a function of its mass M. In turn, requiring 
that MDM makes all the observed DM measured by cosmology, Odm^ 2 = 0.110 ± 0.005 [11], 
we can univocally determine M. As a general rule of thumb, it is well known [1] that ^dm/j 2 « 
3-10~ 27 cm 3 sec -1 / {(JaP) and that a particle with weak couplings a w has a (cta/3) ~ a 2 Mq M that 
matches flj)M for a typical weak scale mass (the so-called WIMP miracle). This is therefore 
what is to be expected for a pure WIMP model such as MDM is. 

More precisely, the computation of the relic abundance has to be performed by solving the 
relevant Boltzmann equation (as we review below) in terms of the detailed annihilation cross 
section of two MDM particles into any SM state. We include the dominant s-wave contribution, 
but also the subdominant p-wave, which yields an 0(5%) correction on the final Qbm and the 
effect of the renormalization of the SM gauge couplings up to the MDM mass scale M (also 
an 0(5%) modification). On top of this standard computation, however, we have to include 
the non-perturbative electroweak Sommerfeld corrections, that have a very relevant effect. In 
fact, this phenomenon significantly enhances non-relativistic annihilations of DM particles with 
mass M > My /a, when they exchange force mediators of mass My with a coupling strenght a. 
In the case of MDM, this is simply the ordinary SM weak force mediated by W ± and Z, and we 
will verify 'a posteriori' that the relation M> M\y±/ot2 is indeed verified for the MDM masses 
M that we will obtain. Therefore the Sommerfeld effect is automatically present in the theory 
and has to be taken into account. We review the basics of the Sommerfeld effect in Appendix 
A. 

Let us now describe in more detail the machinery of the computation and the results in our 
specific case. The generic Boltzmann equation that governs the DM abundance as a function 
of the temperature T reads 

dY /Y 2 \ T r°° / /s\ 

' zhz * i >* ia = «^ L ds s Ki ( T )»"<•> (3) 

where z = M/T, Ki is a Bessel function, Z = (1 — | — ^f)~\ the entropy density of SM particles 
is s = 2tt 2 g* s T 3 / '45, Y = n DM /s where tzdm is the number density of DM particles plus anti- 
particles, and Y e(l is the value that Y would have in thermal equilibrium. The adimensional 
'reduced annihilation cross section' is defined as 

where s, t are the Madelstam variables and the sum runs over all DM components and over all 
the annihilation channels into all SM vectors, fermions and scalars, assuming that SM masses 
are negligibly small. In the case of MDM, we can write a single equation for the total DM 
density, in particular neglecting the small splitting AM that we will discuss in Sec. 3.1, because 
DM scatterings with SM particles maintain thermal equilibrium within and between the single 
components. In this way the formula automatically takes into account all co-annihilations 
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among the multiplet components. We also ignore the Bose-Einstein and Fermi-Dirac factors as 
they are negligible at the temperature T ~ M/26 relevant for DM freeze-out. 

In full generality, for fermionic DM with SU(2)l x U(1)y quantum numbers n and Y we 

get 
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where x = s/M 2 and /3 = y 1 — 4/x is defined in the DM DM center-of-mass frame. The 
first line gives the contribution of annihilation into vectors, the second line contains the sum of 
the contributions of annihilations into SM fermions and vectors respectively. The gauge group 
factors are defined as 
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where the sum is over all SM vectors A = {Y, W 1 , W 2 , W 3 } with gauge coupling generators 
T A . We have defined gx as the number of DM degrees of freedom for a multiplet with Y = 0: 
gx = n for scalar DM, gx = 2n (4n) for fermionic Majorana (Dirac) DM. 

The DM freeze-out abundance is accurately determined by the leading two terms of the ex- 
pansion for small /3, that describe the s-wave and the p-wave contributions. This approximation 
allows to analytically do the thermal average in eq. (3): 
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This general formalism can now be easily reduced to the case of the MDM particle, the 
fermionic 5-plet with Y — 0. The s-wave and p-wave coefficients are simply given by 
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(10) 



The large figures at the numerators simply reflect the 'large' number (n = 5) of components in 
the multiplet: their coannihilations make for an overall large cross section parameter. 

As anticipated, on top of this the annihilation cross sections are enhanced by the non- 
perturbative Sommerfeld effect. The effect is strongly dependent on the velocity /3 (in units of 
c) of the DM particles: the total annihilation cross section aft grows as ft —>■ as illustratated 
in fig. lb. MDM annihilations at the epoch of freeze-out (when ft ~ 0.2) are enhanced by a 
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Figure 1: Left panel: The relic abundance of the MDM 5-plet as a function of its mass M. The 
solid red line corresponds to the full result that includes non-perturbative Sommerfeld enhanced 
annihilations (the dashed red line corresponds to the computation in the SU (2) i- symmetric 
limit), while the dashed blue includes perturbative s-wave and p-wave tree level annihilations. 
The horizontal band corresponds to the measured Udm^ 2 , so that its intersection with the red 
line individuates the DM mass (slightly below 10 TeV) and indicates its uncertainty interval. 
Central panel: Velocity dependence of the Sommerfeld- enhanced MDM annihilation cross 
section. Right panel: The DM DM annihilation cross section in the galactic halo today, 
relevant for indirect detection signatures discussed in Sec. 5. The dashed lines show what the 
perturbative result would be without Sommerfeld corrections. 



factor of a few; astrophysical signals due to MDM annihilations in the galaxy today (where 
(3 ~ 1CT 3 ) will be much more enhanced, up to a few orders of magnitude, as discussed in Sec. 5. 

Solving the Boltzmann equation with all these ingredients allows to compute the relic 
abundance of the MDM particle Odm = Y{z — > oo)sM / p CT - lt as a function of its mass M. 
The numerical result is plotted in fig. la. For a given value of the mass, the impact of non- 
perturbative corrections is relevant, as they enhance the annihilation cross section and there- 
fore reduce the corresponding relic abundance. Matching the relic abundance to the measured 
^dm^ 2 = 0.110 ± 0.005 allows therefore to determine 

M = (9.6 ± 0.2) TeV. (11) 

For comparison, the value without Sommerfeld corrections (blue dashed curve in fig. la) would 
be about 4 TeV. 



In summary: the only free parameter of the theory, the DM mass M, is fixed by matching 
the observed relic DM abundance. Not surprisingly, its value turns out to be broadly in the 
TeV range, because MDM is a pure WIMP model for which the so called 'WIMP miracle' 
applies. The value actually ends up being somewhat higher (eq. (11)) because the 5-plet has 
many components so that coannihilations are important and because Sommerfeld corrections 
enhance the annihilation cross section. 



3.1 Mass splitting 

Supersymmetric models typically feature a model-dependent electroweak-scale mass difference 
between the neutralino DM candidate and its chargino partners: the (typically multi-GeV) mass 
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difference originates through tree-level mass mixings. In the MDM case, instead, at tree level 
all the components of the multiplet have the same mass M computed in the previous section, 
and then one-loop electroweak corrections make the charged components slightly heavier than 
the neutral one (if the contrary were true, one would have a charged lightest stable particle, 
which is of course not phenomeno logically allowed). 

In full generality, for a fermionic or scalar candidate with overall mass M and hypercharge 
Y, the mass difference induced by loops of SM gauge bosons between two components of X 
with electric charges Q and Q' is found to be 
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(12) 
(13) 



^\Jr 2 — 4)/2 and s\y the sine of the weak angle. In the numerically relevant 

r— +0 

limit M 3> Myy t z the one-loop corrections get the universal value f(r) ~ 2nr. For the 
fermionic 5-plet with Y — 0, the splitting between the neutral component (the DM particle) 
and its Q — ±1 partners equals therefore to 



AM = a 2 M w skr 



2 #W 



(166 ± 1) MeV. 



(14) 



In general, the mass splitting between charged and neutral components can be intuitively 
understood in terms of the classical non-abelian Coulomb energy (the energy stored in the 
electroweak electric fields that a point-like charge at rest generates around itself, which can be 
thought as an additional contribution to its mass with respect to the mass of an equivalent 
neutral particle). Indeed, for a scalar or fermion with a gauge coupling g under a vector with 



mass My, the Coulomb energy is: 
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As SXJ(2) L invariance is restored at distances r <^ l/M W)Z , the UV divergent term cancels out 
when computing the correction to the intra-multiplet mass splitting. Therefore we understand 
why the effect in the limit M 3> M\y t z does not depend on the DM spin, and why the neutral 
component is lighter than the charged ones. This intuitive picture allows us to guess that, when 
considering the low-energy DM/nuclei relevant for DM direct detection, the spin-independent 
cross section will be suppressed only by 1/M^ Z (and not by the much smaller 1/M 2 ), because 
what scatters is the cloud of electro- weak electric fields that extend up to r ~ 1/Mw,z- This 
agrees with the Feynman diagram computation presented in the next section. 



4 Direct Detection signatures 

Direct searches for DM aim to detect the recoils of nuclei in a low background detector pro- 
duced by the rare collisions of DM halo particles on such nuclei. So far, the DAMA/Libra 
experiment [13] has reported the detection of an annual modulation of the total number of 
events compatible with the effect that the motion of the Earth in the DM halo would produce. 
The CDMS [9] and Xenon [10] experiments have published exclusion bounds. 



S 




Figure 2: One loop DM / 'quark scattering for fermionic MDM with Y = (two extra graphs 
involving the four particle vertex exist in the case of scalar MDM). 

As discussed in Sec. 2, MDM candidates with Y = have vanishing DMA/" direct detection 
cross sections at tree level (see eq. (17)). The scattering on nuclei J\f proceeds therefore at one- 
loop, via the diagrams in fig. 2 that involve one of the charged components X ± of the multiplets. 
An explicit computation of these one-loop diagrams is needed to understand qualitatively and 
quantitatively the resulting cross section. Non-relativistic MDM/quark interactions of fermionic 
X with mass M 3> My/ 3> m q are described by the effective on-shell Lagrangian 



SIS = (B '-(i ± 2y)')^E [(4 + 4) { * X]m * M - aF^^l^'l 



(15) 



where the + (— ) sign holds for down-type (up-type) quarks q = {u, d, s, c, b, t}, rrih is the Higgs 
mass and m q are the quark masses. The first operator gives dominant spin-independent effects 
and is not suppressed by M; the second operator is suppressed by one power of M and gives 
spin-dependent effects. Parameterizing the nucleonic matrix element as 

{N\Y,m q qQ\N) = fm N (16) 

q 

where is the nucleon mass, the spin-independent DM cross section on a target nucleus M 
with mass M/v is given by 

^(DMV DMAO = („> - ir^MMl (4 + . (17, 

The case of scalar X is not much different: the M-independent contribution to ogi is equal to 
the fermionic result of eq. (17) but there is no spin-dependent effect. 

Assuming = 115 GeV and / m 1/3 (QCD uncertainties induce a one order of magnitude 
indetermination on osi 2 ) we find therefore for the fermionic MDM 5-plet 

o-si = 1.2 ■ 1(T 44 cm 2 . (18) 

As usual [1, 14, 15], cxsi is defined to be the cross section per nucleon. The prediction is a 
definite number (as opposed to the large areas in the plane M/a that is covered by typical 
supersymmetric constuctions by varying the model parameters) and Fig. 3 shows that this 
value is within or very close to the sensitivities of experiments currently under study, such 
as Super-CDMS and Xenon 1-ton [16]. The annual modulation effect of the DAMA/Libra 
experiment [13] cannot be explained by MDM candidates, since they have too large masses and 
too small cross sections with respect to the properties of a WIMP compatible with the effect. 

2 More precisely, one needs to consider the effective Lagrangian for off-shell quarks, finding various operators 
that become equivalent only on-shell. Their nucleon matrix elements can differ; we ignore this issue because 
presently it is within the QCD errors. 
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Figure 3: Parameter space of the spin-independent cross section per nucleon: the 

MDM 5-plet is represented by the green dot. The predictions is univocal (having assumed that 
the matrix element f = 1/3 and 115 GeV as a value for the higgs mass), while typical minimal 
SUSY models span a large portion of the parameter space, suggested here by the area with dotted 
contours. The dashed lines indicate the sensitivity of some future experiments [16]. Other non- 
minimal candidates, discussed in Appendix B, are reported for completeness: the fermionic 
3-plet and the fermion and scalar doublet. 

5 Indirect Detection signatures 

Indirect searches are one of the most promising ways to detect Dark Matter. DM particles in 
the galactic halo are expected to annihilate and produce fluxes of cosmic rays that propagate 
through the galaxy and reach the Earth. Their energy spectra carry important information on 
the nature of the DM particle (mass and primary annihilation channels). Many experiments are 
searching for signatures of DM annihilations in the fluxes of 7 rays, positrons and antiprotons 
and there has been recently a flurry of experimental results in this respect: 

• data from the PAMELA satellite show a steep increase in the energy spectrum of the 
positron fraction e + / (e + +e~) above 10 GeV up to 100 GeV [17], compatibly with previous 
less certain hints from HEAT [18] and AMS-01 [19]; 

• data from the PAMELA also show no excess in the p/p energy spectrum [20] compared 
with the predicted background; 

• the balloon experiments ATIC-2 [21] and PPB-BETS [22] report the presence of a peak 
in the e + + e~ energy spectrum at around 500-800 GeV; 

• the HESS telescope has also reported the measurement [23] of the e + +e~ energy spectrum 
above energies of 600 GeV up to a few TeV: the data points show a steepening of the 
spectrum which is compatible both with the ATIC peak (which cannot however be fully 
tested) and with a power law with index —3.05 ± 0.02 and a cutoff at ~2 TeV. 
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Some words of caution apply to the balloon results: the Monte Carlo simulations that such 
experiments need in order to tag e ± and infer their energy have been tested only up to LEP 
energies; the excess is based on just a few data-points that are not cleanly consistent between 
ATIC-2 and the smaller PPB-BETS; emulsion chambers (EC) balloon experiments [25] do not 
show evidence for an excess, although they have larger uncertainties. The upcoming results of 
the FERMI/LAT mission on the e + + e~ energy spectrum [26] will hopefully soon indicate in 
a more definitive and precise way the shape of the spectrum. 3 

In full generality [27], the PAMELA positron and anti-proton data, if interpreted in terms 
of DM, indicate either (i) a DM particle of any mass (above about 100 GeV) that annihilates 
only into leptons, not producing therefore unseen antiprotons or (ii) a DM particle with a mass 
around or above 10 TeV, that can annihilate into any channel. Adding the balloon data, the 
mass is pinned down to about 1 TeV and only leptonic channels are allowed. 

Minimal Dark Matter has definite predictions for the fluxes of positrons, antiprotons and 
gamma rays, presented in [5] before the announcements of the experimental results (and even 
before some of the relevant experiments, such as FERMI, started taking data). The comparison 
with the data that are now available allows to fix the astrophysical uncertainties and test the 
model. The MDM 5-plet annihilates at tree level into W + W~ , and at loop level into 77, jZ, 
ZZ. The relative cross-sections are significantly affected by the non-perturbative Sommerfeld 
corrections discussed in Sec. 3 and in Appendix A. The best-fit values are: 

3 3 

(crv) ww = 1.1- 10~ 23 — , {<™) 77 = 3- 10~ 25 — (19) 

sec sec 

Annihilation cross sections into 7Z and ZZ are given by 

cr~ /Z = 2cx 77 / tan 2 # w = 6.5<r 77 , a zz = cr 77 / tan 4 # w = 10.8<r 77 (20) 

(that actually hold for any MDM candidate with Y — 0). However, as shown in fig. lc, as a 
consequence of Sommerfeld corrections the DM DM annihilation cross sections exhibit a quite 
steep dependence on Mom and can vary by about one order of magnitude around these central 
values within the range allowed at 3a by the cosmological DM abundance. The cross section 
also depends on the DM velocity v, reaching a maximal value for v — > 0, as shown in fig. lb. 
The average DM velocity in our galaxy, v ~ 10~ 3 , is however low enough that av is close to its 
maximal value, which we assume. 

The energy spectra of e + , p and 7 per annihilation at production, as generated by PYTHIA 
supplemented by appropriate custom routines that allow to take into account the spin correla- 
tions of SM vectors 4 , are represented in fig. 4. We next need to consider where these fluxes of 
particles are produced in the galaxy and how they propagate to the Earth. 

5.1 Positrons 

The positron flux per unit energy from DM annihilations in any point in space and time is 
given by <& e +(t,x, E) = v e +f/4ir (units 1/ GeV ■ cm 2 ■ s • sr) where v e + is the positron velocity 
(essentially equal to c in the regimes of our interest) and the positron number density per unit 
energy, f(t,x,E) = dN e +/dE, obeys the diffusion-loss equation: 

9 J- t - K(E) • V 2 / - A {b{E )f) = Q (21) 



3 Note added. See footnote 5. 
4 More details are given in ref. [5] 
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DM DM -> W+W~ 



DM DM -> ZZ 




Figure 4: Energy spectra of e + ,p, 7 produced by non-relativistic DM DM annihilations into 
W + W~ and into ZZ, the only two relevant channels for MDM (together with 77, which produces 
a 7 line at the MDM mass, and 7Z , which can be deduced from the 77 and ZZ channels). The 
e + have a secondary component (dashed green line shown on the W + W~ plot), that dominates 
at large x ~ 1. 







Positrons 




Antiprotons 
K (kpc 2 /Myr) V COQV (km/s) 




Model 


5 


K (kpc 2 /Myr) 


5 


L (kpc) 


MIN 


0.55 


0.00595 


0.85 


0.0016 13.5 


1 


MED 


0.70 


0.0112 


0.70 


0.0112 12 


4 


MAX 


0.46 


0.0765 


0.46 


0.0765 5 


15 



Table 2: Propagation parameters for charged (anti)particles in the galaxy (from [29], [35]). 



with diffusion coefficient K(E) = K (E/ GeV) 5 and energy loss coefficient b(E) = E 2 /( GeV • 
Te) with te = 10 16 s. They respectively describe transport through the turbulent magnetic 
fields and energy loss due to synchrotron radiation and inverse Compton scattering on CMB 
photons and on infrared galactic starlight. Eq. (21) is solved in a diffusive region with the shape 
of a solid flat cylinder that sandwiches the galactic plane, with height 2L in the z direction 
and radius R = 20 kpc in the r direction [28]. The location of the solar system corresponds 
to x — (tq, Zq) = (8.5 kpc, 0). The boundary conditions impose that the positron density / 
vanishes on the surface of the cylinder, outside of which positrons freely propagate and escape. 
Values of the propagation parameters S, K and L are deduced from a variety of cosmic ray 
data and modelizations. They represent a source of uncertainty over which to scan in order to 
reach the final predictions for the fluxes. We consider the sets presented in Table 2 [29]. Finally, 
the source term due to DM DM annihilations in each point of the halo with DM density p(x) 
is 

where k runs over all the channels with positrons in the final state, with the respective thermal 
averaged cross sections av. Several galactic DM profiles are computed on the basis of numerical 
simulations: isothermal [31], Einasto [32], Navarro- Frenk- White [33], Moore [34] (roughly in 
order of cuspiness at the galactic center). The choice of profile introduces a further element of 
astrophysical uncertainty over which to scan. 
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Figure 5: The astrophysical propagation functions for positrons and antiprotons. 

Left: The 'halo function' I (\d) of eq. (23) that encodes the astrophysics of DM DM annihilations 
into positrons and their propagation up to the Earth. The diffusion length is related to energy 
losses as in eq. (14) of [5], that also provides fit functions for all cases. Right: The somewhat 
analogous p astrophysical function R(T) of eq. (25). In both cases, the dashed (solid) [dotted] 
bands assumes the min (med) [max] propagation configuration of table 2. Each band contains 3 
lines, that correspond to the isothermal (red lower lines), NFW (blue middle lines) and Moore 
(green upper lines) DM density profiles. 



The solution for the positron flux at Earth can be written in a useful semi-analytical form [29, 
30]: 

„ , 1 / n \ 2 /-Mum 

where B > 1 is an overall boost factor discussed below, \d(E,E') is the diffusion length from 
energy E' to energy E. The adimensional 'halo function' I(Xd) [29] fully encodes the galactic 
astrophysics and is independent on the particle physics model. Its possible shapes are plotted 
in fig. 5a for most common choices of DM density profiles and set of positron propagation 
parameters. 

The flux of positrons from DM annihilations has to be summed to the expected astrophysical 
background. We take the latter from the CR simulations of [36] as parameterized in [37] 
by = 4.5£ a7 /(l + 650 E 2 - 3 + 1500 E 4 - 2 ) for positron and = $^,prim + $bkg,scc = 

0.16 E- 1 - 1 / (1 + 11 E - 9 + 3.2 E 2 - 15 ) + 0.70£ a7 /(l + 110 E 1 ^ + 580 E 4 - 2 ) "for electrons, with 
E always in units of GeV. These not-so-recent background computations have recently been 
revised and questioned [38]: background shapes with a downturn around energies of a few 
GeV have been investigated in order to incorporate the PAMELA excess as a feature of the 
background. 

Finally, the DM density in our galaxy might have local clumps that would enhance the 
positron flux by an unknown 'boost factor' B > 1. We take it as energy independent: this is 
a simplifying (but widely used) assumption. Detailed recent studies [39, 40, 41, 42] find that 



13 




Energy in GeV 

Figure 6: Positron fraction. The positron fraction from MDM galactic annihilations is com- 
pared with the data from PAMELA and previous experiments. We have taken a boost factor of 
50 (see text). We show the DM signal (the lines) and the total e + fraction when summed to the 
background (shaded area). The main result (solid line and shaded area) is computed assuming a 
benchmark NFW DM profile and MED propagation parameters. The fainter dashed lines corre- 
spond to changing the propagation parameters to MIN (lower) and MAX (upper). The fainter 
dotted lines correspond to changing the DM profile to isothermal (lower) and Moore (upper). 

a certain energy dependance can be present, subject to the precise choices of the astrophysical 
parameters. Within the uncertainty, these studies also converge towards small values of B 
(except for extreme scenarios), with B = O(10) still allowed. 

On the basis of the ingredients above, the fluxes at Earth of positrons from MDM annihila- 
tions can be compared with the experimental results. This is shown in fig. 6. One sees that the 
predictions agree very well with the PAMELA results on the whole range of energies. We have 
assumed B = 50, which is the value found to provide the best fit to positrons, electrons and 
antiprotons data (discussed below) combined. This is quite a large value, in tension with the 
determinations discussed above. On the other hand, lower values (down to about 20) would 
still give a reasonably good fit and in any case the MDM annihilation cross sections of eq. (19) 
carry a one order-of-magnitude uncertainty. In order to be conservative, we prefer to quote 
the boost value 'as is', instead of looking for possible optimizations. The figure also illustrates 
that the DM signal is only very mildly affected by changing the DM density profile (dotted 
lines). It somewhat depends on the e + propagation model in our galaxy (dashed lines); this 
uncertainty will be reduced by future measurements of cosmic rays and is however present only 
at E < M DM . 

5.2 Electrons + positrons 

The computation of the fluxes of e + + e~ from MDM annihilations is just a rearrangement of the 
calculations for positrons presented in the previous section. Fig. 7 shows the predicted flux as 
compared to the results from the balloon experiments ATIC, PPB-BETS and EC and HESS. It 
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Figure 7: Sum of electrons and positrons. The flux of electrons and positrons from MDM 
annihilations is compared to the data from ATIC, PPB-BETS, EC and previous experiments, 
plus the datapoints from HESS and FERMI. The main result (solid line and shaded area) is 
computed assuming a benchmark NFW DM profile and MED propagation parameters. The 
fainter dashed lines correspond to changing the propagation parameters to MIN (lower) and 
MAX (upper). The fainter dotted lines correspond to changing the DM profile to isothermal 
(lower) and Moore (upper). 

is apparent that the MDM predictions are not compatible with the peak individuated (mainly) 
by the ATIC data points: a spectrum which is flat up to the higher energies, with a smooth 
endpoint somewhat below M would be expected. The HESS datapoints indicate a steepening 
of the spectrum with respect to GeV energies. They are compatible with the shoulder of the 
ATIC peak, which however cannot be fully tested. 

If the presence of the ATIC peak will be confirmed in future data sets, therefore strongly 
indicating a DM mass around 1 TeV, MDM will be falsified. 

Assuming that the astrophysical background is a power-law (HESS data however indicate a 
steepening around 1 TeV), the MDM 5-plet predicts a slightly-harder quasi-power-law e + + e~ 
spectrum up to several TeV energies. 5 

5.3 Antiprotons 

The propagation of anti-protons through the galaxy is described by a diffusion equation anal- 
ogous to the one for positrons. Again, the number density of anti-protons per unit energy 

5 Note added. The ATIC peak, that was incompatible with Minimal Dark Matter, is now contradicted by 
the new more precise FERMI data [66] (unless an exceptionally bad energy resolution is assumed for FERMI), 
that we just superimposed to fig. 7, without modifying the MDM prediction to better fit the new data. The 
problem is now that FERMI data are consistent with the steepening apparent in the HESS data (supplemented 
by the new results at lower energy in [67]). If this feature will be confirmed and cannot be attributed to the 
astrophysical background, Minimal Dark Matter will be excluded as an interpretation of the PAMELA and 
FERMI e ± excesses. 
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p kinetic energy in GeV 

Figure 8: Antiprotons. The antiproton over proton ratio (at top of the atmosphere) from 
MDM annihilations, compared to the recent PAMELA data. We have assumed the same boost 
factor as for positrons (B = 50). The main result (solid line and shaded area) is computed 
assuming a benchmark NFW DM profile and MED propagation parameters. The fainter dashed 
lines correspond to changing the propagation parameters to MIN (lower) and MAX (upper). 
The fainter dotted lines correspond to changing the DM profile to isothermal (lower) and Moore 
(upper). 



f(t, x, T) = dNp/dT vanishes on the surface of the cylinder at z — ±L and r = R. T = E — m p 
is the p kinetic energy, conveniently used instead of the total energy E (a distinction which will 
not be particularly relevant for our purposes as we look at energies much larger than the proton 
mass m p ). Since m p ^> m e we can neglect the energy loss term, and the diffusion equation for 
/is 

df d 

-± - K(T) ■ V 2 / + — (sign(z) / Konv) = Q-2h 8{z) T ann /. (24) 

The pure diffusion term can again be written as K(T) = K Q (3 (p/ GeV) 5 , where p = (T 2 + 
2m p T) 1 / 2 and (3 = Vp/c = (l — m 2 /(T + m p ) 2 ) 1//2 are the antiproton momentum and velocity. 
The V conv term corresponds to a convective wind, assumed to be constant and directed outward 
from the galactic plane, that tends to push away p with energy T < 10 m p . The diferent sets of 
values of the parameters are given in table 2. The last term in eq. (24) describes the annihilations 
of p on interstellar protons in the galactic plane (with a thickness of h = 0.1 kpc L) with 
rate r ann = (n H + 4?^ 3 'nHe)o~pp n Vp, where n H ~ 1/cm 3 is the hydrogen density, n^ e ~ 0.07 n H is 
the Helium density (the factor 4 2//3 accounting for the different geometrical cross section in an 
effective way) and the given esplicitely in [5, 43, 30]. We neglect the effect of "tertiary 
anti-protons" . This refers to primary p after they have undergone non-annihilating interactions 
on the matter in the galactic disk, losing part of their energy. 

In the "no-tertiaries" approximation that we adopt, the solution [44, 45, 46] for the antipro- 
ton flux at the position of the Earth $p(T, r ) = Vp/(4n)f acquires a simple factorized form 
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(see e.g. [35]) 

*^>- s sG^)V)£iM*§ (25) 

x k 

where B is the boost factor. The k index runs over all the annihilation channels with anti- 
protons in the final state, with the respective cross sections; this part contains the particle 
physics input. The function R(T) encodes all the astrophysics and depends on the choice of 
halo profile and propagation parameter set. It is plotted in fig. 5 for several possible choices. 
Finally, for completeness we also take into account the solar modulation effect, due to the 
interactions with the solar wind, that distorts the spectrum via a slight increase of the low 
energy tail, as described in more detail in [47, 5]. 

The astrophysical background is predicted by the detailed analysis in [48], the results of 
which we find to be well reproduced by a fitting function of the form log 10 $^ kg = —1.64 + 
0.07 r - r 2 - 0.02 r 3 + 0.028 r 4 with r = log 10 T/GeV. We take for definiteness the flux 
corresponding to the 'MED' propagation parameters. Particularly favorable is the fact that 
the uncertainty in the estimates of the background is quite narrow around 10 — 100 GeV, where 
results are expected soon. 

Fig. 8 shows the results for the final p flux from MDM annihilations at Earth (at the top 
of the atmosphere), normalized to the one of protons, compared to the background and to 
the most recent PAMELA data (older experimental data point are not anymore significant). 
We have assumed the same boost factor as for positrons, although in principle it could be 
different [40]. As apparent, the agreement with the data is quite good. A prominent excess is 
predicted to show up at energies slightly above those that have been probed by PAMELA so 
far. This is a consequence of the large mass of the MDM candidate and also of the annihilation 
channel into W + W~ . The figure is produced with the benchmark choices of NFW and MED. 
The shape of the spectrum is relatively independent from the propagation model and the halo 
profile. Different p propagation models instead change the overall signal rate by about one 
orders of magnitude. Different halo profiles with fixed p & make only a difference of a factor of 
a few, which can be interpreted in terms of the fact that the signal is not dominated by the far 
galactic center region, where profiles differ the most. 

5.4 Gamma rays 

Gamma ray signals from DM annihilation in the regions where DM is more dense can be very 
significant, but also depend strongly on the astrophysical assumptions, such as in particular the 
galactic DM profile. The spectra of high energy 7 rays from the galactic center in the Minimal 
DM case have been presented in [4, 5]: the spectrum is characterized by a (smeared) line at E w 
M and a continuum that extends to lower energies. Most astrophysical detections of gamma 
ray fluxes are however compatible with power law fluxes: the most conservative approach is 
therefore to assume that the observed fluxes are of standard (yet unknown) astrophysical origin 
and impose that the flux from DM does not excede such observations. This imposes stringent 
constraints on the DM annihilation cross section and the DM halo profiles. 

A complete, model independent calculation has been carried out in this respect in [49], 
which considers bounds from high energy gamma rays from the galactic center, the galactic 
ridge and the Sagittarius dwarf galaxy, but also radio waves from the galactic center (produced 
by synchrotron radiation in the strong magnetic field by the electrons and positrons from DM 
annihilations). Such bounds apply of course in particular to the case of a particle with a 9.6 
TeV mass and W + W~ main annihilation channel (fig.4 in [49]) such as Minimal DM: it is 
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found that the annihilation cross section in eq. (19) is (marginally) allowed in the case of a 
NFW DM profile in the Milky Way and the Sagittarius dwarf galaxy. It has however to be 
assumed that no boost factor is present for signals from the galactic center, which is possible 
if DM clumps are tidally disrupted in the central regions of the Milky Way. Alternatively, 
perhaps more realistically, but in tension with state-of-the-art numerical simulations, choosing 
less steep profiles such as isothermal allows to pass the constraints. Recent studies [50, 51] 
along these lines for dwarf galaxies find comparable or looser constraints. 



6 Collider searches 

At an accelerator like the Large Hadron Collider (LHC), which will soon(er or later) operate 
colliding pp at y/s = 14 TeV, DM particles can in principle be pair produced, and this is a very 
promising way to search for DM and measure its properties. 

In the case of MDM, for any multiplet with Y = and arbitrary n the partonic total cross 
sections (averaged over initial colors and spins) for producing any of its component are 

« _ 9 - _ 9 ~ _ 9x9t(n 2 - 1) j f3 2 if X is a scalar . . 

<?ud - °du - ^uu ~ - 13g24 ^ P ■ | 3 _ ^2 if ^ is a f ermion 

where the subscripts denote the colliding partons, and (3 = \J\ — AM 2 / s is the DM velocity with 
respect to the partonic center of mass frame. Production of non-relativistic scalars is p-wave 
suppressed in the usual way. It is obvious that the production of the 5-plet with M = 9.6 TeV 
is kinematically impossible. Possible upgrades of LHC luminosity and magnets are discussed 
in [52]. 

If produced in collisions, MDM has a clean signature: the small mass splitting AM among 
the DM components makes charged MDM component (s) enough long-lived that they manifest 
in the detector as charged tracks. Irrespectively of the DM spin the life-time of DM 1 * 1 particles 
with Y = and n = {3,5, 7, . . .} is r ~ 44cm/(n 2 — 1) and the decay channels are precisely 
determined as 



DM ± 


-> DMV 


: T n = 


(n 2 - 1) 
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-> DM°e ±( i7 ) e 


■ r e = 
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: r M = 


o.i2 r e 



f3 /2 / 777,^ 



4tt V AM 2 



BR n = 97.7% 



Gp AM 5 BR e — 2.05% (2?) 



60tH 



BR M = 0.25% 



having used the normalization f n = 131 MeV [53] and the AM of eq. (14), which accidentally 
happens to be the value that maximizes BR„-. The DM + life-time is long enough that decays 
can happen inside the detector. On the contrary, the faster decays of DM ±=t particles (present 
for n > 5) mostly happen within the non- instrumented region with few cm size around the 
collision region. Measurements of r and of the energy of secondary soft pions, electrons and 
muons constitute tests of the model, as these observables negligibly depend on the DM mass M. 
On the contrary, measurements of the total number of events or of the DM velocity distribution 
would allow to infer its mass M and its spin. Although SM backgrounds do not fully mimic 
the well defined MDM signal, at an hadron collider (such as the LHC) their rate is so high that 
it seems impossible to trigger on the MDM signal. 

Notice that extra SU(2)l multiplets that couple (almost) only through gauge interactions 
tend to give a LHC phenomenology similar to the one discussed above, irrespectively of their 
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possible relevance for the MDM problem. On the contrary, DM candidates like neutralinos are 
often dominantly produced through gluino decays, such that DM is accompanied by energetic 
jets rather than by charged tracks. 

7 Other phenomenological signatures 

The characteristic small mass splitting AM between the neutral DM and the charged compo- 
nents of the multiplet allows for a few other interesting albeit quite speculative manifestations 
of MDM. 

7.1 Accumulator 

One could envision accelerating a large amount of protons p or nuclei TV in an accumulator ring 
and having them collide on the DM particles of the galactic halo bath, that would therefore 
act as a diffuse target. At energies E PiA r > AM the CC collision DMA/" — > A/' ± DM =F becomes 
kinematically possible and the DM T product would scatter out of the beam pipe giving a 
signature (collider experiments need instead a much larger energy E P >M in order to pair 
produce DM). An estimate of the event rate gives 

dN_ _ Pdm _ 10 N p p DM TeV a 
dt~ £ v ° M ~ £ yearl0 20 0.3GeV/cm 3 M 3cr 1 ' 

where pdm is the local DM density and e is the detection efficiency, related e.g. to the fraction 
of beam that can be monitored. <jq = GpM^/ir = 1.1 10 _34 cm 2 is a reference partonic cross 
section: at large energies Myy <E P <M the MDM-quark cross section approaches 3<7o(n 2 — 1)/4. 
Proton drivers currently planned for neutrino beam experiments can produce more than 10 16 
protons per second, and accumulating N p ~ 10 20 protons is considered as possible, so that 
the number of expected events looks not unreasonable. The main problem however seems 
disentangling the signature from the beam-related backgrounds, such as collisions on residual 
gas. 

7.2 Ultra High Energy Cosmic Rays 

If MDM particles are somehow a component of the Ultra High Energy (UHE) Cosmic Rays, 
the already relatively long half-life of DM 1 * 1 is Lorentz-dilated to macroscopical distances, and 
the possibility opens of seeing spectacular signatures in the form of km-long charged tracks of 
DM 1 * 1 in detectors like Icecube or Antares. Consider a flux of DM particles that is crossing the 
Earth at UHE. Via Charged Current (CC) interactions with nucleons of Earth's matter, DM 1 * 1 
particles are produced. Being charged particles traveling in a medium, these loose a part of 
their energy and eventually decay back to DM particles. This chain of production and decay is 
analogous to the process that tau neutrinos undergo in matter ("^-regeneration"). In the n = 5 
case DM ±=t states are also produced and decay. Such a system is described by a set of coupled 
integro-differential equations for the evolution of the fluxes of DM , of DM ± and of DM ±=t as 
a function of their energy and the position through the Earth. At high energies, the CC cross 
section interactions are large (for E > 10 15 eV DM particles interact on average at least one 
time along the whole diameter of the Earth) , but so are also the energy losses and the Lorentz- 
dilated mean life of DM 1 * 1 and DM ±± , so that a detailed analysis (numerical or semi-analytical) 
is necessary. The results are reported in fig. 9b: the relevant phenomenological variable is the 
fraction of particles, at a given energy, that on average would emerge from their Earth crossing 
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Figure 9: MDM in UHECR. The fraction of DM particles that exit from an Earth-crossing 
journey in the charged states DM* 1 and DM ±ziz (dashed lines; the results of numerical and semi- 
analytical computations are represented by the paired lines that almost coincide), as a function 
of the exit energy E out . The red dot-dashed line show the ratio E out /E- in and allows to reconstruct 
the initial energy with which they must have entered the other side of the Earth. For instance, 
one reads that about 10% of the DM particles that exit the Earth at 10 19 eV do so in the DM ± 
state; they have lost about 2/3 of their E m in the journey. This fraction has to be multiplied by 
the flux of UHECR and by the fraction of it which is assumed to be made of MDM particles in 
order to get an actual event rate. UHECR have been detected up to few -10 20 eV. 



in the charged state and therefore would leave detectable tracks in the neutrino telescopes. 
One sees that, at high energies, such a fraction can be sizable. However, at those very high 
energies the overall flux of Cosmic Rays is not vanishing but fairly faint, so that the overall 
number of events would be limited. Moreover, it seems challenging to construct a mechanism 
for accelerating MDM particles at UHECR energies in the first place. They could however be 
a component of cosmic rays generated by the decay of ultra-heavy particles, in the so called 
'top-down' scenario for the explanation of the origin of UHECR. Finally, particle identification 
in a detector like IceCUBE or Antares poses difficulties: DM 1 * 1 would roughly look like muons 
(produced by neutrino interactions) with fake energy = E±Pcc/P^ ~ 10~ 4 -E±, because muon 
energy losses are approximatively given by —dE^fdt = a + (3^E^ with (3^ ~ 0.2/kmwe [59] {(3 
is roughly inversely proportional to the particle mass). One therefore needs to carefully study 
the energy loss profile along the charged tracks to see the difference between a DM ± with 
E± ~ 10 18 eV and a muon with E^ ~ 10 14 eV. Larger and more densely instrumented detectors 
seem necessary to study these issues. 



8 Conclusions and outlook 

As experimental searches for Dark Matter proceed more vigorously than ever, many models may 
finally have the chance of being tested. Supersymmetric DM, Kaluza-Klein DM, Little Higgs 
DM etc. all make predictions that however (i) inscribe themselves into the context of oddly not- 
yet-found EW-scale New Physics; (ii) typically depend on a large number of unknown model 
parameters, so that DM phenomenology remains dark; (iii) usually rely on the existence of 
some extra feature (R-, KK-, T-parity) imposed by hand in order to be kept stable against 
decay into SM particles. 

The Minimal Dark Matter proposal takes on a more... minimal approach: focussing on the 
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Table 3: Main phenomenological properties of the preferred MDM candidate. Estimates of the 
theory uncertainties on the numerical values are given in the text. 

DM problem only, we add to the Standard Model just one extra multiplet X with electroweak 
interactions (such that it does not ruin the successful predictions that the SM makes in its 
renormalizable limit: conservation of baryon number, lepton number, lepton flavor, etc.), and 
investigate whether it can constitute a good DM candidate: electrically neutral, stable, pro- 
duced in the right amount via the thermal freeze-out in the Early Universe and not excluded 
by direct DM detection results. We find that indeed the construction selects one preferred 
candidate, the fermionic SU(2)l quintuplet with hypercharge Y — 0, the main phenomenolog- 
ical properties of which are listed in Table 3. These properties are univocally computed, as 
no free parameters are present in the model. Also, no extra feature is introduced by hand to 
guarantee its stability: simply no decay modes exist consistently with the SM gauge symmetry 
(analogously to the stability of the proton in the SM). Phenomenological predictions can be 
univocally listed: 

— The univocal MDM signature at colliders is production of DM 1 * 1 , that manifests as a non- 
relativistic charged track (straight despite the magnetic field B ~ Tesla in the detector) 
that decays with a relatively long life-time of r ~ 1.8 cm (due to the small splitting) 
into DM°7r ± , leaving a quasi-relativistic curved track. However, this signature will not be 
visible at LHC, because the mass of the preferred MDM candidate is too large (see [52] 
for possible LHC upgrades) and because this signature is too complex for triggers (given 
their speed and the rate of QCD backgrounds). 

+ The next generation of direct detection experiments, such as Super-CDMS or Xenon 
1-ton, prospect to be sensitive to the spin-independent scattering cross section of the 
MDM candidate on nuclei, see fig. 3. MDM cannot account for the controversial annual 
modulation claimed by DAMA/Libra. 

+ Indirect DM searches constitute arguably the most interesting testing ground for the 
MDM model. The predictions of the theory had been presented in [5] and can now be 
confronted to the recent experimental results from PAMELA, ATIC and HESS, of course 
under the assumption that these are interpreted in terms of DM annihilations. The main 
features of the predicted spectra are determined by the following main properties of the 
model: 

• the DM mass is very large (9.6 TeV), so that fluxes are expected to extend to multi- 
TeV energies; 

• the predominant annihilation channel is into W + W~ channel, that produces fluxes 
of all species (e + , p, 7 and u) with a characteristic spectrum; 
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• the total annihilation cross section is very large (of the order of 1CT 23 cm 3 /sec) 
thanks to the non-perturbative Sommereld enhancement which is present due to the 
exchange of EW gauge bosons between the annihilating particles, as reviewed in 
Appendix A. 

The profile of DM distribution in the galaxy and the propagation of the fluxes of charged 
anti-matter in the galaxy introduce an uncertainty on the final spectra that is however of 
limited impact, essentially because the large mass of the MDM particle produces fluxes 
of high energy final products that do not travel for long in the galaxy otherwise they 
would loose all their energy (the dependence on the propagation parameters remains 
more significant, at the level of one order of magnitude, on p fluxes). Fig.s 6, 7 and 8 
illustrate that: 

+ The MDM 5-plet prediction agrees with the PAMELA data on the positron 

fraction, if an astrophysical boost factor of B ~ 50 is introduced (reduced to B ~ 5 
within the 3a allowed uncertainty range for the annihilation cross section). The 
model predicts a continuous rise up to an energy corresponding to the DM mass of 
about 10 TeV. 

+ The MDM 5-plet prediction agrees with the PAMELA data on antiproton/proton 
fluxes, assuming the same astrophysical boost factor. The reason of the agreement 
is that the excess in p corresponding to the excess in e + starts at energies slightly 
higher than those probed by PAMELA so far. Therefore the model predicts a very 
relevant anomaly in future PAMELA or AMS-02 data with respect to the expected 
background. 

— The MDM 5-plet flux of the sum of electrons and positrons is not compatible with 
the peak suggested by ATIC and PPB-BETS data at around 600 GeV: an e + + 
e~ spectrum which is flat up to larger energies is expected instead. The HESS 
observations hint to a steepening in the e + + e~ spectrum, without however testing 
the peak. New data should soon clarify the experimental situation, possibly falsifying 
MDM. 6 

Moreover: 

~ The production of gamma rays (and of radio synchrotron emission from the e ± ) from 
the MDM 5-plet at the galactic center and in dwarf galaxies are compatible with the 
observations for an NFW or Einasto profile if the astrophysical boost is assumed not 
to be present for observations in the gamma channel towards the galactic center, or 
if a not-too-steep DM profile such as isothermal is assumed (somewhat in tension 
with numerical simulations). 

~ Neutrinos from the annihilations of MDM trapped at the Sun's center are not ex- 
pected to be detectable, as the fluxes depend on the inverse of the square of the 
DM mass (which is large) and the spin dependent capture cross section on solar 
nuclei is small. Neutrinos from the galactic center can be relevant for detectors with 
sensitivity larger than the current one by a few orders of magnitude. 

6 Note added. The ATIC peak, incompatible with MDM predictions, is now contradicted by the new 
FERMI data [66] (unless an exceptionally bad energy resolution is assumed for FERMI). However FERMI data 
are consistent with the steepening apparent in the HESS data (recently supplemented by the new results at 
lower energy in [67]). 
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~ Anti-deuterium d is formed when the DM annihilation products contain a p and a 
n with momentum difference (in their rest frame) less than about 160 MeV [54]. 
Ref. [55] computed the d spectrum from MDM in the standard spherical approxima- 
tion [54] , finding a yield well below the background. However Minimal Dark Matter 
annihilates into W boosted by a Lorentz factor 7 = M/Mw m 120, and the jet 
structure of the products has to be taken into account. This has been done in [56]. 
The total d yield is the same as the one from decays of a W at rest, but everything is 
boosted by a factor 7, resulting (for the central value of av we assume) into a d flux 
of about 10~ 10 per m 2 ,secsrGeV at energies of about 100 GeV per nucleon. This 
should be then comparable to or above the expected astrophysical d background. 
These large energies are however challenging for the purposes of detection. 

If the anomalies in the fluxes of leptons (e + fraction and e + + e~) in the ranges of energies 
currently explored will turn out to be of astrophysical origin, e.g. from one or more 
local pulsars [57, 58], the predictions of MDM (as well as of most DM models) will be 
drowned in an unsuppressible background. The MDM model is peculiar in predicting, 
however, excesses that continue up to 10 TeV in these channels. Also, the prediction for 
the flux of anti-protons (much more difficult to mimic with astrophysics) would stay, and 
it cannot be reduced below a minimum level: with unit boost factor, the MDM flux is 
still marginally above the background for most choices of parameters, albeit concentrated 
at very high energies. 

+ Constraints on the amount of energy deposited by DM annihilations at the time of BBN 
(when the velocity of the DM particles can be < 10~ 3 ) impose an upper bound on the 
annihilation cross section. They were recently reconsidered by [60] after [61]: the annihi- 
lation cross section of MDM eq. (19) is found to be allowed (see e.g. fig. 2 of [60]; note 
that for MDM av is already at saturation at (3 ~ 10~ 3 , see fig. lb). 

Further experimental results are expected soon. 



A Sommerfeld effects 

The DM DM annihilation cross section gets enhanced if DM particles have a non-relativistic 
velocity (3 <C 1 and there is a long-range attractive force between them. Technically, the 
effect arises because annihilating DM particles cannot be approximated as plane waves, so 
their wave-function must be computed by solving the Schroedinger equation in presence of the 
long-range potential V. For a single abelian massless vector with potential V = a/r the result 
is a = SVperturbative where the Sommerfeld correction is [62] 

S(x) = x = -. 29 

v ; 1 - e nx (3 y ' 

Here a < describes an attractive potential that leads to an enhancement S > 1, and a > 
describes a repulsive potential that leads to S < 1. The above results holds for the s-wave 
(L = 0) partial wave, that dominates in the non-relativistic limit: the s-wave annihilation 
cross section a/3 roughly grows as 1/(3 for (3<a. Higher waves (L > 0) are enhanced by higher 
powers of a/ (3 [64] such that, as (3 — > all partial cross sections grow as 1/(3, but with negligible 
coefficients suppressed by a 2L . 

The Sommerfeld enhancement is automatically present in the MDM case, since the MDM 
particle is coupled to the SM W, Z vectors and it is significantly heavier than them. As a 
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consequence the s-wave annihilation cross section a[3 grows as 1/(3 for 1 > (3 > Mw/M ~ 10 3 . 
This is why a possibly large positron excess had been predicted in the MDM framework before 
the PAMELA results. 7 

In the MDM framework the force vectors (the weak gauge bosons) are massive and non- 
abelian: the computation of the Sommerfeld effect is more involved as one must classify the 
two-body DM DM states according to their conserved quantum numbers: angular momentum 
S, total spin L and total electric charge Q. Ref. [4] presents all the details of the calcula- 
tions in the different cases. Here we focus on the DM°DM° state of the 5-plet, relevant for 
astrophysical signals, which lies in the Q = 0, L = 0, S = sector containing 3 states: 
{DM ++ DM , DM + DM~, DM°DM }. The 'potential' V and the 'annihilation rate' T become 
3x3 matrices, whose off-diagonal components (and especially their signs) do not have an intu- 
itive meaning and are defined in terms of the real and imaginary part of the two-body matrix 
propagators. The explicit result is [4] 

++/8A-4A -2B 
V= + -2B 2A-A 
\ -3V2B 

where A = 166 MeV, A = a em /r + a2C^ !V e~ Mzr jr and B = a2e~ Mwr jr. The Sommerfeld 
enhancement can now be computed by numerically solving the matricial Schroedinger equation. 

Co-annihilation with all other DM components enter in the computation of the cosmological 
freeze-out abundance, so that a lengthy computation is needed. We here discuss how one can 
perform a simplified computation in the limit of unbroken SU(2)l, neglecting the SM vector 
masses M^^ and the MDM intra-multiplet mass splitting A. This approximation is good 
enough for the cosmological computation. Group theory allows to reduce the unbroken non- 
abelian Sommerfeld enhancement to a combination of abelian-like enhancements in the various 
sectors [63], classified according to the conserved quantum numbers: L,S and iso-spin /. We 
can focus on s-wave annihilations (L = 0) and on states with / < 5, that annihilate at tree level 
into two SM particles. The Pauli principle restricts the allowed states to be (/, 5*) = (3, 1), (1, 0) 
and (5, 0). After performing the relevant thermal average, the s-wave coefficient of eq. (10) gets 
Sommerfeld-enhanced to 

c s = ^S{-Sa 2 ^Mjf) + ^<S(-3a 2 yM/T) + (1 + h^S(-Qa 2 y/M/T) . (31) 
^jr ^ , ^ , v 24 tt ^ , 

(I,S) = (l,0)->W a W a (I,S)=(5,0)->W a W b (/,s)=(3,i)^*« , ,ra 

The resulting DM freeze-out abundance is shown in the left panel of fig. 1 (dashed line), 
compared with the full computation (continuous line). 
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B Next-to-minimal Minimal Dark Matter 

In this Appendix we briefly sketch the directions into which the model can be extended, if the 
request of full minimality is abandoned and some tools of DM model building (commonly used 
elsewhere) are introduced. 

7 In other DM models, ad hoc light new particles are often introduced such that the Sommerfeld enhancement 
allows to accomodate the PAMELA anomaly (as it suggests a a (5 at (3 ~ 10 -3 which is a few orders of magnitude 
larger than the af3 w 3 10 26 cm 3 /sec at f3 ~ 0.2 suggested by the cosmological DM abundance). 
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In Sec. 2 we rejected the candidates which have dimension-4 or dimension-5 interactions with 
SM fields and can therefore decay into SM particles in a time much shorter than the age of the 
universe. By invoking some ad hoc extra symmetry that prevents the decays, these candidates 
can of course be reallowed. This is actually also what happens in the case of the best know 
WIMP DM candidates: in SUSY models, by assigning odd i?-parity to the supersymmetric 
partners including the LSP DM, the decay operators, that feature just one DM state, are 
forbidden. An equivalent mechanism is used in little-Higgs models (T-parity), in 'universal' 
extra dimension models (KK-parity), etc. Some of the MDM candidates which are rescued 
by the imposition of the extra symmetry resemble to particles that appear in a variety of 
other contexts: e.g. scalar triplets in little-Higgs models; fermion or scalar triplet in see-saw 
models; KK excitations of lepton doublets or of higgses in extra dimensional models; higgsinos, 
sneutrinos, winos in supersymmetric models. In general, however, these models possess a 
plethora of other parameters so that typically their DM candidates behave very differently 
from the MDM candidates. But in some corner of the parameter space the phenomenological 
features can coincide. The signatures of some of these candidates (in particular the 'wino-like' 
fermionic triplet and the scalar triplet) have been studied in [4, 5]. 

Speaking of decays, we note that the main MDM candidate (the fermion 5-plet, cosmologi- 
cally stable in the minimal construction) acquires a non-negligible decay mode if A is somewhat 
below the Planck scale. The effective dimension-6 operator that dominantly contributes to the 
decay is XLHHH* / A 2 . It induces 4-body decays, such as DM -> FW£Z L Z L (where L 
denotes longitudinal polarization), with rate T ~ M 5 /4-7rA 4 . Such decays can provide an alter- 
native interpretation of the PAMELA excess: A w 10 16 GeV gives a life-time of about 10 25 sec. 
2-body decays, such as DM — > £ ± H^ =F have a rate T 2 ~ v 4 M/ (4-7r) 5 A 4 which is subdominant as 
long as¥> 4irv. 

As discussed in Sec. 4, MDM candidates with 7^0 have an elastic cross section which is 
2 3 orders of magnitude above present direct detection bounds [9] due to the exchange of 
a Z boson and have therefore been rejected. A possible way to reinclude them is to abandon 
minimality by introducing another state that mixes with the DM: this has the effect of splitting 
the components of the Dirac DM particle by an amount 8m such that the lightest one becomes 
a Majorana fermion which cannot have a vector-like coupling to the Z boson. This is actually 
what happens in the well-known case of SuSy Higgsino DM, via the mixing with Majorana 
gauginos. In our case, if a Sm larger than the DM kinetic energy and smaller than Mq — M is 
generated, the only consequence would be to suppress the direct detection DM signals and the 
rest of the phenomenology would be essentially unchanged. 

In the case of scalar DM, additional operators not listed in eq. (1) are generically present: 
4-scalar interactions of X with higgses and X quartic self-interactions: 

& D -c\ H {X*T x X) (H*T£H) - c\' H \X\ 2 \H\ 2 - — {X*T x Xf - ^\X\\ (32) 

where Tjj are SU(2)l generators in the representation to which R belongs and c = 1 (1/2) for a 
complex (real) scalar. In the strictly minimal theory these terms have been assumed vanishing, 
as they introduce the new unknown couplings A#, X' H , \x and X' x . But if minimality is relaxed 
they can be introduced and they allow to enlarge the parameter space of the theory. These extra 
interactions do not induce DM decay (because two X are involved, and we assume (X) = 0). 
They affect instead the computations of the mass splitting, of the direct detection inelastic 
cross section and of the relic density (therefore the determination of the DM mass). Indeed, 
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the coupling Xh splits the masses of the components of X by an amount 



AM ^|Arl| iw 

AM M v 1 

having inserted (H) = (v , 0) with v = 174 GeV and AT| = 1. This cannot be neglected with 
respect to the splitting which is separately generated by loop corrections (discussed in Sec. 3.1) 
as soon as Xh ~ 0.01. Also, Xh and X' H open the possibility of extra annihilations into higgses 

|A' g | 2 + (n 2 -l)|A g | 2 /16 
((T ^ )cxtra = 16. M 2 ^ (34) 

(no interference terms are present). The contribution from Xh is relevant for Xh ~ 0.01, the 
same value that also produce a non-negligible mass splitting. The contribution from \' H is 
relevant for |A^| ~ gy-,g\ or larger. In these regimes, to compute how much the inferred value 
of M would be affected one simply adds the contribution in eq. (34) to those in eq. (6): in 
general, the increase of the total annihilation cross section implies the increase of the inferred 
MDM mass, and the mass of the low n candidates are more affected than the high n ones, 
because the gauge-mediated annihilation cross sections are relatively less important for the 
formers. For example, a large A'^ = 1 would increase the predicted value of M by a factor 
2.4 for n = 2, by 20% for n = 3, by 2% for n = 5, and by 0.5% for n = 7. Note that 
Xj{ ~ g\ and X' H ~ g\ are the values predicted by some solutions to the hierarchy problem, 
such as supersymmetry and gauge/higgs unification. Finally, non-zero couplings Xh and X' H 
can also produce an elastic cross section cxsi on nuclei, much larger than the one in eq.(4), as an 
UV-divergent effect that corresponds to a renormalization of |rY| 2 |if | 2 operators is generated. 

We note that when X is a neutral scalar singlet, the non-minimal annihilations in eq. (32) are 
the only existing ones. In this case, the observed amount of DM is obtained for M « 2.2 TeV| A'^ | 
(we are assuming M 3> M z ; for generic values of M the correlation between M and X' H was 
studied in [65]). 

Finally, one can envision a scenario in which more than one MDM multiplet is present at 
the same time. Or also a scenario in which several distinct copies ('flavors', in analogy with 
the SM families) of the same MDM multiplet exist. In both cases, to a good approximation 
their abundances evolve independently, such that the observed DM abundance is reproduced 
for lower values of the respective DM masses. In this more general situation, the M values in 
table 1 must be reinterpreted as upper bounds on M. This might bring some of the candidates 
into the reach of the LHC. 
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